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t-h Abstract 

The pairing gaps, heat capacities and level densities are calculated within the BCS-based quasi- 
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,_L| particle approach including the effect of thermal fluctuations on the pairing field within the pairing 
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model plus noncollective rotation along the z axis for 60 Ni and Ge nuclei. The analysis of the 
numerical results obtained shows that, in addition to the pairing gap, the heat capacity can also 



J> serve as a good observable to detect the appearance of the pairing reentrance in hot rotating nuclei, 
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whereas such signature in the level density is rather weak. 
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I. INTRODUCTION 



In the collective rotation of a deformed nucleus, the rotation axis is perpendicular to the 
symmetry axis, and the Coriolis force, which breaks the Cooper pairs, increases with the 
total angular momentum so that at a certain critical angular momentum all Cooper pairs 
are broken. The nucleus undergoes then a phase transition from the superfluid phase to the 
normal one (the SN phase transition). This is called the Mottelson-Valatin effect fTJ. A 
similar effect is expected in spherical nuclei, although the rotation is no longer collective. 
The total angular momentum is made up by those of the nucleons from the broken pairs, 
which occupy the single-particle levels around the Fermi surface and block them against 
the scattered pairs. The pairing correlations decrease until a sufficiently large total angular 
momentum M c , where the pairing gap A completely vanishes. At finite temperature (T 7^ 
0), the increase of T relaxes the tight packing of quasiparticles around the Fermi surface, 
which is caused by a large angular momentum M > M c , and spreads them farther away from 
the Fermi level. This makes some levels become partially unoccupied, therefore, available 
for scattered pairs. As the result, when T increases up to some critical value Ti, the pairing 
correlations are energetically favored, and the pairing gap reappears. As T goes higher, the 
increase of a large number of quasiparticles eventually breaks down the pairing gap at T 2 
(> Ti). This phenomenon, predicted by Kamuri (2j and Moretto [3j, is called thermally 
assisted pairing or anomalous pairing, and later as pairing reentrance by Balian, Flocard 
and Veneroni jl]. 

However, it has been shown already in the 1960s that the sharp SN-phase transition at 
M = M c in the Mottelson-Valatin effect is an artifact of the BCS method. As a matter of 
fact, a proper particle-number projection before variation has removed the discontinuity in 
the pairing gap as a decreasing function of the angular momentum [5] . Similarly, by taking 
the effect of thermal fluctuations in the pairing field into account, the SN phase transition 
predicted by the BCS theory is smoothed out. The gap A(T) of a non- rotating nucleus does 
not collapse at T c ~ 0.568A(T = 0), but monotonically decreases with increasing T, remain- 
ing finite even at T ^> T c j6j [7] . This result is reconfirmed by shell-model calculations of 
pairing energy as a function of excitation energy j8], and by embedding the exact eigenvalues 
of the pairing problem into the canonical ensemble j9]. By considering an exactly solvable 
cranked deformed shell model Hamiltonian it has also been shown that the pairing gap, 
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quenched at T = and high rotational frequency, reappears at T\ (<C T c ) at M > M c [TO] . 
However, different from the prediction by the BCS theory, the pairing gap does not vanish 
at T > T v 

The behavior of hot rotating nuclei can be put in correspondence with superconductors in 
the presence of an external magnetic field, where the magnetic field plays the role as that of 
the nuclear rotation. The reentrance of superconducting correlations, which is also known 
as the unconventional superconductivity, has been the subject of recent theoretical and 
experimental studies in condensed matter. The Grenoble High Magnetic Field Laboratory 
has recently discovered that URhGe becomes superconducting at low temperature in the 
presence of a strong magnetic field (between 8 and 13 T), well above the field value of 2 
T, at which superconductivity is first destroyed [UJ. The reentrance of superconductivity 
under the magnetic field is interpreted as to be caused by spin reorientation [12] , which bears 
some similarity with the reappearance of the scattered pairs in rotating hot nuclei discussed 
above. 

Recently, we have developed an approach based on the finite-temperature BCS (FTBCS) 
that includes the effects due to quasiparticle-number fluctuations in the pairing field and 
the z projection of angular momentum at T ^ 0, which we call as the FTBCS1 (with "1" 
denoting the effect due to quasiparticle-number fluctuations) [14]. This approach reproduces 
well the effect of smoothing out the SN phase transition at T ^ as well as the pairing 
reentrance in hot (noncollectively) rotating nuclei. In the latter case, for M > M c the 
pairing gap also reappears at T\ and remains finite at T > T\. It has also been pointed 
out that the microscopic mechanism of the nonvanishing gap at high T is the quasiparticle- 
number fluctuations, which are ignored in the conventional BCS theory. A refined version 
of the FTBCS 1 also includes the contribution of coupling to pair vibrations within the 
self-consistent quasiparticle random-phase approximation [14] . 

That the pairing reentrance is not an artifact of the mean field analysis, but a robust 
physical effect, has been obtained by exact diagonalization of the 2D attractive Hubbard 
model, where a nonmonotonic filed dependence of the pair susceptibility in the presence of 
the external magnetic field was found for various cluster sizes both in the weak and strong 
coupling limit [13] . Nonetheless, the experimental extraction of the pairing gap in hot nuclei 
is not simple because one has to properly exclude the admixture with the contribution of 
uncorrelated single-particle configurations from the odd-even mass difference [15]. Therefore 
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the detection of the pairing reentrance effect by using the experimentally extracted pairing 
gaps seems to be elusive, especially when the formally derived pairing gap has a value smaller 
than the average spacing between the single-particle levels. 



Meanwhile, the heat capacity has been extracted from the experimental level densi- 
ties [16J. The existence of a bump or an S shape on the curve of the heat capacity at 
T ~ T c allows one to discuss about the smoothing of the SN phase transition in finite nuclei. 
In a recent calculation of the heat capacity in 72 Ge within the Shell Model Monte Carlo 
(SMMC) approach, by reconfirming the pairing reentrance effect, the authors of Ref. [17] 
claimed that they found a local dip in the heat capacity at rotation frequency of 0.5 MeV 
at T ~ 0.45 MeV, and a corresponding local maximum on the temperature dependence of 
the logarithm of level density. They associated such irregularities in the heat capacity and 
level density as the signatures of the pairing reentrance. There are, however, two concerns 
regarding these results. The first one is that, as well-known, at such low temperature the 
SMMC approach produces quite large error bars. As a consequence, instead of approaching 
zero as it should be to fulfill the third law of thermodynamics, the SMMC heat capacity at 
T < 0.5 MeV jumps to 25 ~ 30 (Fig. 4 of [17J), which makes the statement on the signature 
of the pairing reentrance ambiguous. The second one is that the SMMC in Ref. [TT] used 
the same Fock-space single-particle energies of the shells (0/lp — 0gld2s) for both neutron 
and proton spectra. Hence, the difference between neutron and proton spectra came solely 
from the difference in the valence particle numbers outside the closed-shell core of 40 Ca, 
which ignored altogether the Coulomb barrier in the proton spectrum. Our FTBCS1 theory 
is free from such deficiency at low T, and it works well with the schematic as well as single- 
particle spectra, which are obtained from the realistic Woods-Saxon potential. Therefore, 
in the present paper we will calculate the heat capacity as well as the level density within 
the FTBCS1 theory to see if these quantities can be used to identify the pairing reentrance 
phenomenon in realistic nuclei at finite temperature and angular momentum. 

The paper is organized as follows. The formalism for the calculations of thermodynamic 
quantities such as pairing gap, heat capacity, and level density in hot non-collectively rotating 
nuclei within the FTBCS and FTBCS1 theories is presented in Sec. [TT| The results of 



numerical calculations are analyzed in Sec. Ill The paper is summarized in the last section, 
where conclusions are drawn. 
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II. FORMALISM 



We consider the pairing Hamiltonian describing a spherical system rotating about the 
symmetry z axis [3J: 

H = Hp — XN — 7M , (1) 

where A and 7 are the chemical potential and rotation frequency, respectively. Hp is the 
standard pairing Hamiltonian of a system, which consists N particles interacting via a 
monopole pairing force with the constant parameter G (the BCS pairing Hamiltonian), 
namely 

Hp = ^ e k{a\ k a+k + aL k a -k) ~ a W-k a -k' a k' , (2) 

k kk' 

with a± k (a±k) being the creation (annihilation) operators of a particle (neutron or proton) 
with angular momentum k, projection ±m k , and energy e k . The particle-number operator 
N and total angular momentum M, which coincides with its z projection, are given as 

N = ) j (a\_ k a +k + a ] _ k a^ k ) , M = ' S ^m k {a ] +k a +k - a ] _ k a^ k ) . (3) 

k k 

After the Bogoliubov transformation from the particle operators, a\ and a k , to the quasi- 
particle ones, a\ and a k , 

a\ = u k a\ + v k a_ k , a_ fc = u k a_ k - v k a\ , (4) 

the Hamiltonian ([T| is transformed into the quasiparticle one "H, whose explicit form can be 
found, e.g., in Refs. [HHH]. 

As has been discussed in Refs. [21 [3J [H], for a spherically symmetric system, the 
laboratory-frame z axis, which is taken as the axis of quantization, can always be made 
coincide with the body-fixed one, which is aligned with the direction of the total angular 
momentum within the quantum mechanical uncertainty. Therefore the total angular mo- 
mentum is completely determined by its ^-projection M alone. For systems of an axially 
symmetric oblate shape rotating about the symmetry axis, which in this case is the principal 
body-fixed one, this noncollective motion is known as "single-particle" rotation. The pairing 
reentrance effect was originally obtained within the BCS theory in Refs. [21 E] by considering 
such systems described by Hamiltonian Its physical interpretation based on the thermal 
effect, which relaxes the tight packing of quasiparticles around the Fermi surface due to a 
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large angle momentum M > M c , and spreads them farther away from the Fermi level, fits 
well in the framework of this "single-particle" rotation. However, as has been pointed out 
in Ref. [3], for non-spherical nuclei, and specifically in the case of axially symmetric ones, 
the spin and angular momentum projections on the symmetry axis are not good quantum 
numbers. In this case the formalism used here is not completed because it does not include 
the angular momentum's component perpendicular to the symmetry axis. Cranking model 
might serve as a better solution of the problem in this situation. This remains to be investi- 
gated because the results obtained within the Lipkin model with J x cranking did not reveal 
any pairing reentrance so far [2T] . On the other hand, in the region of high level densities 
(at high excitation energies and/or high T) the values of angular momentum projection on 
the symmetry axis will be mixed among the levels, which worsen the axial symmetry. The 
melting of shell structure will also eventually drive nuclei to their average spherical shape. 

A. FTBCS1 equations at finite angular momentum 

The FTBCS1 includes a set of FTBCS-based equations, corrected by the effects of 
quasiparticle-number fluctuations, for the level-dependent pairing gap A*, average particle 
number N, and average angular momentum M. The derivation of the FTBCS1 equations 
was reported in detail in Ref. [133, so we d° present here only the final equations. The FT- 
BCS1 equation for the pairing gap is written as a sum of two parts, the level-independent 
part A and the level-dependent part SA^, namely 



A fc = A + 6A k , 



(5) 



where 




G 




—u k v k , 



(6) 



where 




(7) 



with the quasiparticle-number fluctuations 5N"l at nonzero angular momentum 




(8) 



6 



The corrections due to coupling to pair vibration beyond the quasiparticle mean field at 
finite temperature and angular momentum are significant only in light nuclei like oxygen 
or neon isotopes, whereas they are negligible for medium and heavy nuclei (See Figs. 6 - 
8 of Ref. |14j). As the present paper considers two medium nuclei, 60 Ni and 72 Ge, these 
corrections on the FTBCS1 equations are neglected in the numerical calculations. 

The equations for the particle number and total angular momentum are the same as Eqs. 
(25) and (26) of Ref. [14J, namely 



N 



2 £ 



1/ t 



M 



E 



m AH -n fc . 



(9) 



The system of coupled equations ^ - ^ are called the FTBCS1 equations at finite angular 

momentum. Once the FTBCS1 equations are solved, the total energy £, heat capacity C 

and entropy S of the system are calculated 

d£_ 
df ' 



£ = (H) , C 



S = ~ I>* lim fc + " < M 1 - < ) + n k ln % + (1 - n k ) ln (! - n k )] • (10) 



B. Level density 



Within the conventional FTBCS, the level density is calculated as the invert Laplace 
transformation of the grand partition function [3] 



p(E, N, M) = j^ry J d/3 j da <j> dpe 



(11) 



where a = /3X, \i = (3^, and S is entropy of the system. The saddle-point approximation 



gives a good evaluation of the integral ( 11 ). As the result, the total level density of a system 



with N neutrons and Z protons is given as 



where S = Sn + Sz and 



D 



p(E,N,M) = 


e s 




(2vr)V 


D ' 


d 2 n 


d 2 a 








dapjdaz 


da^fdfi 




d 2 n 


a 2 n 


d 2 n 


d 2 n 


da.zda.fq 


da% 


dazdfj, 


dazdp 




d 2 n 


d 2 n 


d 2 n 


dfidajsf 


dfidaz 


dfj. 2 


8^8/3 


d 2 n 


d 2 n 


d 2 n 


d 2 Q 


dfidaff 


dpdaz 




d/3 2 



(12) 



(13) 
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The logarithm of the grand-partition function of the systems is given as 



D, = tt N + tt z = S + a N N + a z Z + fj,M — f3£ . (14) 

The derivation of Q with respect to a and fi can be seen explicitly in Eqs. (25)-(35) of Ref. 
[3]. Within the FTBCS1, the grand-partition function has the same form as that given by 



Eq. (14) of the FTBCS. Therefore, the first and second derivatives of the FTBCS1 grand- 
partition function are the same as those of the FTBCS ones. The only difference comes 
from the first derivatives of the pairing gap with respect to a, //, and /3 because of the 
quasiparticle- number fluctuations in the FTBCS 1 gap equation ([6]). In this case, instead of 
the simple Eqs. [(33)-(35)] of Ref. [3], the derivatives become rather complicate expressions, 
which are obtained by taking the first derivatives of the left and right-hand sides of Eq. ^ 
with respect to a, fi, and (3. Amongst the three derivatives of the FTBCS1 gap, dA k /df3 
can be obtained by using its definition, namely 

dA k _ 2 dA k _ ^ 2 A fc (T + £T)-A fc (r) 

dp ~ 1 dT ~ 1 6T ' (i5j 

and can be easily calculated numerically by choosing an appropriate value of 5T as the 
input parameter. The other two derivatives, dA k /da and dA k /dfi, must be calculated from 
their explicit analytic expressions because a and (5 are two Lagrange multipliers, which are 
obtained by solving the FTBCS 1 equations. The final equations for dA k /da and dA k /d[i 
are derived as 

E (rf-t + B *) + ( c ? + °r - §) f'-t + w + ^ • < 16 > 

k ' N ' 

E (^^f + B> ) + + D ' k ~ §) 4^ + E ' k = ' (17) 

where 

A k = 4 {(A - e fc ) 2 (l - n+ - n k ) + PA 2 k E 2 k [nt + n k - 2(n+) 2 - 2{n k f]} , (18) 
^k 

Bk = - ^^^ {l-ni-n k - PEKnt + n k - 2«) 2 - 2(n k f) , (19) 

Ct = n ^Sl^lps K A - <*) 2 (1 % ) " ^lE k [(nt - I) 2 + n k (2n+ - n k )]} , 

1 1 n k n k ) ^k 



(20) 



C: 



[l-nt-n- k YEl 



d: 



d; 



{(A - e k ) 2 (l -n+- n k ) - 0AlE k [{n k - l) 2 + n + k (2n k - n+)]} , 

(21) 

nt ~ n; + PE k {(ni - l) 2 + n k (2nt - n k ))} , (22) 



+ (l-n+)A fc (A-e fc ) 
{l-n+-n k YEl 1 

*n ~"P A "-i~J k) {l -nt- n k + PE k [{n k - I) 2 + ni(2n k - „+)]} , (23) 



A' 



- -3 {(A - e fc ) 2 (l -n+- n k ) + /3A 2 E fc [n+(l - n+) + n» (1 - nj" )]} 
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, (3m k A h 



CI 



Ek 
(A - e k f 



it 



+(l-nt)+n k (l-n k )} 



^-<-n k )El 



+(l-nt)+n k (l-n k )] 



[l-nt-n k YEl 



K(l - nt)[{nt - l) 2 + n~(2n+ - n~ k )\ 



+ n~ k (1 - n fe )[(rifc - l) 2 + n+(2n fe - n+)]} 
K = -^^[nt(l-ni) + n k (l-n k )] 



(24) 
(25) 
(26) 



(27) 
(28) 



By solving first the FTBCS1 equations, then Eqs. (16) and (17), one obtains dA k /da and 
dA k /da as functions of T at a given value of the total angular momentum M. 



III. ANALYSIS OF NUMERICAL RESULTS 

The numerical calculations are carried out for two realistic 60 Ni and 72 Ge nuclei. The 
latter is considered in order to have a comparison with the results obtained within the 
SMMC approach in Ref. [IT]. The single-particle spectra for these two nuclei are obtained 
within the axially deformed Woods-Saxon potential [22], whose parameters are chosen to be 
the same as those given in Ref. [18]. All the bound (negative energy) single-particle states 
are used in the calculations. The quadrupole deformation parameters (3 2 are equal to and 
-0.224 for 60 Ni and 72 Ge, respectively. The pairing interaction parameters are adjusted so 
that the pairing gaps at T = fit the experimental values obtained from the odd-even mass 
differences. These values are Gn = 0.347 MeV, which gives A^v = 1.7 MeV for neutrons in 
60 Ni and Gn = 0.291 MeV, which gives A at = 1.7 MeV for neutrons in 72 Ge. For protons 
in 72 Ge, the value Gz = 0.34 MeV is chosen to give Az = 1.5 MeV, whereas there is no 
pairing gap for the closed-shell protons (Z = 28) in 60 Ni. Because the pairing gap ^ is 
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level-dependent, the level-weighted gap A is considered, which is defined as A = A^/fi 
with the total number Q of levels in the deformed basis [In the case of spherical basis it 
becomes A = £\(2j + 1)A,/ ^.(2j + 1) [14j]. 

Shown in Fig. [I] are the neutron level-weighted pairing gap A, the heat capacity C, and 
the ratio C /T obtained as functions of T at several values of the total angular momentum 
M. The left column represents the predictions by the standard FTBCS, whereas the results 
obtained within the FTBCS1 are displayed in the right column. Both approaches show the 
pairing reentrance in the gaps at M = 4 and 6 h, namely the gap increases with T up to 
T ~ 0.3 MeV, then decreases as T increases further. Because of the quasiparticle-number 
fluctuations, the FTBCS1 gap does not collapse at T c as the FTBCS one, but decreases 
monotonically at high T. At M = 14 h, while the FTBCS gap completely vanishes at all T, 
the FTBCS1 gap shows a spectacular reentrance effect, namely it increases from the zero 
value at T = up to around 0.3 MeV at T ~ 0.8 MeV, and then slowly decreases as T 
further increases. 

The heat capacities obtained within the FTBCS and FTBCS1 look alike, except for the 
region around T c , where the quasiparticle-number fluctuations smooth out the sharp SN 
phase transition so that the sharp local maximum is depleted to a broad bump. In the 
region, where the pairing reentrance takes place, namely at T ~ 0.3 MeV and M = 4 or 
6 h, a weak local minimum is seen on the curve representing the temperature dependence 
of the heat capacity similarly to the feature reported in Ref. [T7j. This local minimum is 
magnified by using the ratio C/T so that the latter might be useful in experiments as a 
quantity to identify the pairing reentrance. However, when the gap is too small as in the 
pairing reentrance at M = 14 h, the heat capacity C {C/T) obtained within FTBCS1 is 
almost identical to that predicted by the FTBCS, where the gap is zero. 

For 72 Ge, both neutron and proton gaps exist, which cause two peaks in the temperature 
dependence of the heat capacity obtained within the FTBCS, as shown in Fig. [2] (c). 
The overall features of C and C/T for 72 Ge are similar to those obtained for 60 Ni. As 
compared with the results of Ref. [T7], where the same single-particle energies in the (0/lp — 
0gld2s) shells were used for both neutrons and protons, and where the pairing reentrance 
was predicted for neutrons, no pairing reentrance effect for neutrons is seen in the results of 
our calculations. On the other hand, the pairing reentrance takes place for protons at M > 
6 h, as shown in Fig. [2] (b) and [2] (f). 
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FIG. 1. (Color online) Level-weighted neutron pairing gap A [(a), (d)], heat capacity C [(b), (e)], 
and heat capacity divided by temperature C/T [(c), (f)] for 60 Ni obtained at different values of 
angular momentum M as functions of T. Panels (a) - (c) show the FTBCS results, whereas the 
predictions by the FTBCS1 are displayed in (d) - (f). 
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FIG. 2. (Color online) Level-weighted pairing gaps for neutrons [(a), (e)], protons [(b), (f)], heat 
capacity C [(c), (g)], and heat capacity divided by temperature C/T [(d), (h)] for 72 Ge obtained 
at different values of angular momentum M as functions of T. Panels (a) - (d) show the FTBCS 
results, whereas the predictions by the FTBCS1 are displayed in (e) - (h). 
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FIG. 3. (Color online) Level densities as functions of T at several values of total angular momentum 
M obtained for 60 Ni [(a), (b)] and 72 Ge [(c), (d)] within the FTBCS (left panels) and FTBCS1 
(right panels). 



Another experimentally measurable quantity, which may help to identify the pairing 
reentrance effect, is the level density. In fact, the authors of Ref. [T7] claimed that the 
pairing reentrance causes an irregularity in a shape of a small local maximum at low T 
on the curve, which describes the temperature dependence of the level density. The level 
densities obtained at several values of the total angular momentum M for 60 Ni and 72 Ge 
are displayed in Fig. [3] as functions of T. These results show a trend of transition of the 
level density from a convex function of T to a concave function after the pairing reentrance 
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FIG. 4. (Color online) Level averaged neutron and proton gaps Ajv,z, heat capacity C, C/T, 
excitation energy E* and level density for 72 Ge as functions of T at several values of M [shown in 
(a)] obtained within the FTBCS1 in the test calculations by using the shells (0/lp — 0gld2s) atop 
the Ca core with the values of Woods-Saxon neutron single-particle energies adopted for both 
neutron and proton spectra. 

occurs. This is particularly clear for 60 Ni by comparing the FTBCS1 predictions for the 
level density at M < 14 h, which are convex functions of T, with that obtained at M — 
14 h, which is a concave function of T. For 72 Ge this trend is less obvious because of the 
existence of proton and neutron pairing gaps with different values of T c within the FTBCS. 
However, contrary to the result shown in the inset of Fig. 4 in Ref. p2], no pronounced local 
maximum that might correspond to the pairing reentrance is seen here in the temperature 
dependence of the level density for 72 Ge. Since the results for the pairing gap, heat capacity, 
and level density strongly depend on the selected single-particle energies, the irregularity 
seen in the temperature dependence of the level density at oo — 0.5 MeV in the inset of Fig. 
4 of Ref. [T7] might well be an artifact caused by using the same single-particle energies for 
both neutron and proton spectra. 

To show that it is indeed the case, we carried out the test calculations by using only the 
(0/lp— 0gld2s) shells on top of the 40 Ca core with the values of Woods-Saxon neutron single- 
particle energies adopted for both neutron and proton spectra of 72 Ge, as in Ref. [IT]. The 
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difference now solely comes from that between the numbers of valence neutrons and protons 
(20 valence neutrons and 12 valence protons). The results of these test calculations within 
the FTBCS1 are shown in Fig. |4j They show the pairing reentrance in the neutron pairing 
gap instead of the proton one, at low M. This is in qualitative agreement with the pairing 
reentrance predicted for neutrons in Fig. 2 of Ref. [IT], where no proton pairing reentrance 
is seen up to u = 0.5 MeV. In our calculations, however, the proton pairing reentrance takes 
place at rather high M = 22 h. In either case, the pairing reentrance is so strong that causes 
the excitation energy E* to decrease slightly with increasing T at low T. This violates the 
second law of thermodynamics. As a consequence, at M > 6 h, the heat capacity becomes 
negative at T < 0.4 MeV. The results for the level density remain essentially the same as 
compared to those previously obtained by using all proton and neutron single-particle levels, 
but with only bound-state single particle energies [Fig. [3] (d)], and no irregularities such 
as a pronounced local maximum are found. In our opinion, the pick on the dotted curve 
in the inset of Fig. 4 in Ref. [T7] emerges because of the two lower values of hip at T ~ 
0.42 and 0.45 MeV. These two lower values are the results obtained by calculating the level 
density in the canonical ensemble p(E) = (3e s /\/2irC, making use of the two large values 
of C equal to around 16 and 28 (with large error bars). However, these large values of the 
heat capacity at low T are the artifacts of the SMMC calculations because the heat capacity 
must be zero (or very small) at T = (or very low T) to avoid an infinite (or very large) 
entropy, which would violate the third law of thermodynamics. Therefore, we conclude that 
the neutron pairing reentrance effect, reported in Ref. [T7J, is caused by the use of the same 
single-particle spectrum for both protons and neutrons, whereas the irregularity seen on the 
curve of lnp in the inset of Fig. 4 of Ref. [UJ is caused by unphysically large values of the 
heat capacity at low T in the SMMC technique. 

The present calculations within the FTBCS and FTBCS1 do not take into account the 
effects of residual interactions beyond the monopole pairing one. It is well known that 
these effects are responsible for strong collective motion in finite nuclei, which leads to the 
increase of nuclear level density. In spherical nuclei the collective enhancement of level 
density is caused by vibrational excitations, whereas in deformed nuclei it comes from the 
collective rotation. The contribution of collective motion to the increase of nuclear level 
density has been studied in detail by Ignatyuk and collaborators starting from the early 
1970s [231 121] • Because the Hamiltonian used in the SMMC calculations of Ref. [T7] included 
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the quadrupole-quadrupole interaction, let us estimate the effect of the collective quadrupole 
vibration on the increase of level density. By using the adiabatic approximation (3) of Ref. 
[2~4] for the enhancement coefficient K vib due to quadrupole vibration, and the experimental 
energies E(2±) of the lowest quadrupole excitation in Ref. [25], we found K v n, ~ 1.06 
and 1.94 for 60 Ni at T = 0.3 MeV and 72 Ge at T = 0.4 MeV, respectively. These values 
of T are those, at which the pairing reentrance starts to show up in these nuclei. With 
the deformation parameter /3 2 = -0.224 adopted in the present calculations for 72 Ge and 
by using Eq. (9) of Ref. [26J, we found the enhancement coefficient K co u{(32) — 1.96 for 
72 Ge at T = 0.4 MeV, whereas for the spherical nucleus, 60 Ni, K co u = K vib = 1.06 at 
T = 0.3 MeV. Therefore, for both nuclei, 60 Ni and 72 Ge, one can expect that the collective 
quadrupole enhancement of level density is not dramatic at the value of temperature, where 
the pairing reentrance is supposed to take place. The contribution of collective motion 
generated by higher multipolarities to the increase of level density is expected to be much 
smaller. In Ref. [2Z] the quasiparticle Tamm-Dancoff Approximation, which includes the 
isoscalar quadrupole-quadrupole interaction and J x cranking, was used to calculate the level 
density within the microcanonical ensemble. The authors of Ref. (27] found 4 < K rot < 6 
and 1.002 < K mb < 1.012 at excitation energy 3 < E < 8 MeV (i.e. at around 0.38 < T < 
0.63 MeV) for 162 Dy. They also found a monotonic decrease of the average pairing gap with 
increasing the excitation energy up to E = 8.5 MeV (T ~ 0.65 MeV), i.e. much higher than 
T c ~ 0.34 MeV. These results are in good qualitative agreement with our estimations. 

Finally, it is worth noticing that the inclusion of the approximate particle-number projec- 
tion within the Lipkin-Nogami method does not significantly alter the behavior of the paring 
reentrance obtained within the FTBCS1 theory (See Fig. 6 of Ref. [H]). This does not 
diminish the value of an approach based on exact particle-number and angular momentum 
projections. In Ref. [28] the exact solution of the nuclear shell model is used to study the 
SN phase transition including residual interactions other than the pairing one. The results 
of Ref. [2H], which fully respect the particle number and angular momentum conservations, 
confirm the presence of a long tail of pair correlations far beyond the BCS phase transition 
region in agreement with the prediction by the FTBCS1. The approach of Ref. [28] does 
not use any external heat bath, which determines the temperature of thermal equilibrium. 
Therefore the nuclear temperature can only be extracted from the level density by using the 
Clausius definition of thermodynamic entropy. This task is not easy because of the discrete 
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and finite nuclear spectra (See, e.g. Ref. [15] and references therein). Nonetheless, instead 
of using temperature, it would be interesting to see if the pairing reentrance takes place in 
the pair correlator as a function of excitation energy at various values of angular momentum 
within the method of Ref. |28|. 



IV. CONCLUSIONS 



The present paper studies the temperature dependences of the heat capacity and level 
density in hot medium-mass nuclei, which undergo a noncollective rotation about the sym- 
metry axis. The numerical calculations, carried out by using the realistic Woods-Saxon 
single-particle energies for 60 Ni and 72 Ge within the FTBCS and FTBCS1 theories, have 
shown the pairing reentrance in the pairing gap at finite angular momentum M and tem- 
perature T. Instead of decreasing with increasing T, the gap first increases with T then 
decreases at higher T. It is demonstrated that the heat capacity C, or rather C/T, and level 
density p can be used to experimentally identify the pairing reentrance effect. The pairing 
reentrance, when it occurs, leads to a clear depletion in the temperature dependence of the 
heat capacity, whereas the level density weakly changes from a convex function of T to a 
concave one. 

Regarding the appearance of the local minimum in the heat capacity because of the 
pairing reentrance, the results of the present paper agree with that of the SMMC calculations 
in Ref. [17]. However, the present results show no pronounced local maximum in the 
temperature dependence of the level density. The pairing reentrance is seen in the proton 
pairing gap of 72 Ge at low M, whereas Ref. [T7J reported this effect in the neutron pairing 
energy. The test calculations by using the same single-particle configuration as that used 
in Ref. [T7J, but obtained within the Woods-Saxon potential, reveals that the neutron 
pairing reentrance in 72 Ge is an artifact, which is caused by the use of the same single- 
particle spectrum for both protons and neutrons, whereas the irregularity on the curve for 
the logarithm of level density, reported in Ref. [T7J , is caused by unphysically large values 
of the heat capacity at low T in the SMMC approach. 
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